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Gravitation governs the expansion and fate of the universe, and the growth of large 
scale structure within it, but has not been tested in detail on these cosmic scales. 
The observed acceleration of the expansion may provide signs of gravitational laws 
beyond general relativity. Since the form of any such extension is not clear, from 
either theory or data, we adopt a model independent approach to parametrising de- 
viations to the Einstein framework. We explore the phase space dynamics of two key 
post-GR functions and derive a classification scheme and an absolute criterion on 
accuracy necessary for distinguishing classes of gravity models. Future surveys will 
be able to constrain the post-GR functions' amplitudes and forms to the required 
precision, and hence reveal new aspects of gravitation. 
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1. Introduction 

Gravitation is the force that dominates the universe, from setting the overall ex- 
pansion rate to forming the large scale structures of matter. Yet our first precision 
tests of gravitation on large scales indicate our understanding of gravity acting on 
the components we see and expect - baryonic and dark matter - is insufficient. The 
cosmic expansion is not decelerating as it should given these ingredients, but accel- 
erating, pointing to either an exotic component with negative active gravitational 
mass (the sum of the energy density plus three times the pressure) or new aspects of 
gravity. The growth and clustering of galaxies likewise do not agree with a universe 
possessing only gravitationally attractive matter within general relativity. 

Having been surprised in our first two tests of cosmic gravity, we naturally look 
to explore expansions of the classical framework. Canonical general relativity (GR) 
with an exotic dark energy component, whether Einstein's cosmological constant 
or otherwise, is certainly one possibility. Extending our knowledge of gravity is 
another, and is what this article concentrates on. 

One question is how to systematically analyse extensions to the known frame- 
work. This can be done by working within an alternate, fully formed theory, but such 
first principles theories are scarce to nil (but see other articles within these Transac- 
tions for some possible guiding principles). Instead we take a phenomenological, or 
(ideally) model independent, approach in analogy to the manner in which compact 
source gravity uses the parameterised post- Newtonian formalism. Such model inde- 
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Functions 



Parameterisation 



Reference 



g, V 

M = 25 - V, ^ = 
PPF: fa = g-^-l,g = 

MGCAMB: 7 : 

s = M = V 



2S 
V 



23 



Bins in fc, z 

yL,u^ ^ a" or bins in z 

model dependent 

E = 1 + Esfl' 
PGA 



A* = 1 + IJ-sO. 



Daniel & Linder 2010 

See translation table in Daniel et al. 2010 

Hu & Sawicki 2007, Hu 2008 

Zhao et al. 2009, Bertschinger & Zukin 2008 

Song et al. 2010, Song et al 2010b 
Zhao et al. 2010 



Table 1. Translation between several different parameterisations of extended gravity and 
the light / growth functions Q and V. 

pendence may put us in good stead to catch signs revealing the underlying nature 
of surprising observations. 

Building on a robust framework of the conservation and continuity equations, we 
parameterise functions from the equations of motion most closely tied to photon 
and matter density perturbation observablcs. This is equivalent to starting from 
the metric potentials themselves, although not directly from an action. In Sec. [5] we 
describe the parameterisation, and consistency structure of the system of equations, 
in more detail. Discussing the types and reach of data enabled by future surveys such 
as BigBOSS, Euclid, and WFIRST, in Sec. [3] we explore the degree of constraints 
that may be placed on these extended gravity, or post-GR, quantities in comparison 
to our knowledge today. We draw connections between this model independent 
approach and representatives of strong coupling and dimensional reduction classes 
of gravity (/(i?) and DGP, respectively) in Sec. SI as well as exploring a different 
tack to distinguishing between gravity models through a phase space analysis similar 
to that used for dark energy. 



2. Framing Gravity 

In the equations of motion for cosmological perturbations in a homogeneous, isotropic 
background four quantities enter: the time-time and space-space metric potentials 
(equal to each other within GR), the mass density perturbation field, and the veloc- 
ity perturbation field (taking a perfect fluid, e.g. ignoring pressure and anisotropic 
stress). Conservation of stress-energy gives the continuity equation relating the den- 
sity and velocity fields, and the Eulcr equation relating the velocity and time-time 
metric potential, so we are left with two free connecting equations. These can be 
chosen, for example, to be the gravitational slip between the two metric potentials 
and a modified Poisson equation relating the space-space potential to the density 
field, or they could be two Poisson-like equations relating one potential to the den- 
sity field and the sum of the potentials to the density. The latter choice turns out to 
give greater complementarity between the parameterised functions, with one closely 
tied to the growth of matter structures and the other closely related to photon per- 
turbations such as gravitational lensing defiection and the integrated Sachs- Wolfe 
(ISW) effect. 

Tabic [1] shows several different forms of model independent parameterisation, 
and how they translate one into the other, adapted and extended from Daniel et 
al. 2010. 

Recently, several groups (Song et al. 2010, Daniel et al. 2010, Zhao et al. 2010, 
Daniel & Linder 2010) have advocated the "light/growth" forms we will use here. 
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due to their close relations with observables and their near orthogonality. The defin- 
ing equations are 



where is the timc-tinic metric potential, (f) is the space-space metric potential 
(in conformal Newtonian gauge), Ga? is Newton's constant, pm is the homoge- 
neous matter density. Am the perturbed matter density (gauge invariant), k is the 
wavenumber, and a is the scale factor. The functions Q{k,a) and V(fc, a) gcner- 
ically are length scale and time dependent. Arising from the sum of potentials, 
Q (meant to evoke an effective Newton's constant) predominantly governs photon 
perturbations through light deflection and the ISW effect. Coming from the veloc- 
ity equation, V predominantly governs growth and motion of structure. They thus 
probe reasonably distinct areas of extension to standard gravity. 

These functions also map onto cosmological observations in somewhat orthogo- 
nal ways. Cosmic microwave background (CMB) data is mostly sensitive to the ISW 
effect and hence Q (since we want the gravitational modifications to be responsible 
for current acceleration, we assume at high redshift, e.g. at CMB last scattering, 
the theory acts like general relativity without modifications). Weak gravitational 
lensing involves Q through the light deflection law, and V to some extent through 
the growth of structure. Galaxy distributions and motions are most sensitive to V. 
Note that the phenomenological gravitational growth index parameter 7 of Linder 
2005, Linder & Cahn 2007 is directly related to V (Daniel et al. 2010). 

Cosmological data, now or in the near future, will not have sufficient leverage to 
reconstruct the general functions G{k, a) and V{k, a). Just as with the dark energy 
equation of state w{a) describing the cosmic expansion, one can only constrain a 
very limited number of parameters describing the functions. Some assume a par- 
ticular time dependence, and possibly neglect scale dependence (e.g. Song et al. 
2010b). More generally, one can use principal component analysis to determine the 
best constrained eigenmodes of the functions (Zhao et al. 2009b). Here we will use 
a similar but simpler approach of dividing the functions into bins of redshift and 
wavenumber, since this provides a more direct physical interpretation of the results: 
gravity is modified in a certain way at low/high redshift and smaller/larger scales. 
We find that two bins in redshift and two in wavenumber, for each of the two post- 
GR functions (we call this the 2x2x2 gravity model), is the extent of the leverage 
that next generation surveys will provide, so more complex parameterisations are 
not useful. 



To effectively constrain the post-GR parameters of G{ki,Zi) and V{ki,Zi), where 
i ~ 1,2 represent the two bins, we need observational data that is sensitive to both 
the effects on the photons (for Q) and the matter growth (for V). Distance measures 
such as Type la supernova distances or baryon acoustic oscillation scales arc useful 
for determining background quantities such as the matter density Qm that might 
have covariance with the post-GR parameters. 

For the current state of the art data we can consider CMB photon perturbation 
spectra (WMAP: Jarosik et al. 2010), supernova distances (AmanuUah et al. 2010), 



k^{(l) + tp) 
—k^Tp 



SttGnO^ PniAjn X Q 
A-KGNO^PmAm X V . 



(2.1) 

(2.2) 
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galaxy clustering (Reid et al. 2010), weak gravitational lensing (CFHTLS: Fu et al. 
2008, COSMOS: Massey et al. 2007), and CMB temperature-galaxy crosscorrelation 
(Ho et al. 2008, Hirata et al. 2008). The results, discussed in detail in Daniel et al. 
2010, Daniel & Linder 2010 (also see Bean & Tangmatitham 2010, and Lombriser 
et al. 2009, Thomas, Abdalla, & Weller 2009 for DGP constraints, Lombriser et al. 
2010 for f{R) constraints), show that while Q is currently bounded to lie within 
10-20% of the GR value for each of the four combinations of low/high wavenumber 
and low/high redshift, V is only weakly limited to within ±1 of the GR value. This 
indicates that current growth probes do not have sufficient leverage to look for 
deviations from general relativity on cosmic scales. Moreover, the two weak lensing 
data sets do not agree with each other, with COSMOS showing consistency with 
GR while CFHTLS gives up to 99% cl deviations at high wavenumbers and low 
redshift. This may be due to difficulties in extracting accurate weak lensing shear 
measurements on small angular scales where the density field is more nonlinear. 

To place significant limits on V and growth, future data sets including galaxy 
clustering and weak lensing surveys covering much more sky area, accurately, and 
ideally to greater depth are required. (Another, nearer term probe will be measure- 
ment of the CMB lensing deflection field.) Surveys such as BigBOSS (Schlegel et 
al. 2009, Stril, Cahn, & Linder 2010), Euclid (Refregier et al. 2010, Martinclh et 
al. 2010), KDUST (Zhao et al. 2010b), LSST (LSST 2009), and WFIRST (Gehrels 
2010) offer great potential gains. Clear understanding of redshift space distortions 
would enable probes of the matter velocity field in addition to the density field and 
serve as a method to measure the growth rate (see below), giving further windows 
on extensions to gravitation theory. 

Figure [1] illustrates one example of future leverage possible on the 2x2x2 
post-GR model-independent parameterisation testing gravitation on cosmic scales. 
The constraints on the growth and V tighten to the few-10% precision level, giving 
tests of 8 different post-GR variables (and all their crosscorrelations) to better than 
10%. That could deliver strong guidance on the nature of cosmic gravity. 

To compactify all the information on testing gravity into a single variable (which 
is not always desirable), we can also examine the gravitational growth index param- 
eterisation approach and the constraints that future data will be able to place on 7. 
Remember that this is only a partial characterisation of extensions to gravity, but 
can serve as an alert to deviations from general relativity (or to matter coupling) 
if the derived value of 7 shows time or scale dependence or is inconsistent with 
"fGR ~ 0.55. 

One promising method to measure 7 is redshift space distortions in the galaxy 
power spectrum (Linder 2008, Guzzo et al. 2008). This depends on the growth rate 
/ = dhiA„i/dlna ~ fi„i(a)''', as well as the growth itself A,„(a) and the galaxy 
bias b{a). One needs strong knowledge of the bias and modeling of the redshift 
distortion form (beyond linear theory) , or excellent data (clear angular dependence 
maps or higher order correlations) to separate out 7 without assuming a form for 
the bias. If the shape of the bias is fixed, keeping the amplitude as a fit parameter, 
then next generation galaxy surveys such as BigBOSS, Euclid, or WFIRST can 
measure 7 to ^ 7%, simultaneously with fitting the expansion history and neutrino 
mass effects on growth (Stril, Cahn, & Linder 2010). 

Another probe is weak gravitational lensing, either by itself (in which case the 
mass is measured and galaxy bias does not enter) or in crosscorrelation with the 
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Figure 1. Filled contours show 68% and 95% cl constraints on V — 1 and Q — 1 for the 
two redshift and two wavenumber bins using mock future BigBOSS, Planck, and WFIRST 
supernova data. The dotted contours recreate the 95% cl contours from Figs. 8 of Daniel & 
Linder 2010 using current data, to show the expected improvement in constraints. The x's 
denote the fiducial GR values (note the offset of current contours may be from systematics 
within the CFHTLS weak lensing data). Adapted from Daniel & Linder 2010. 



galaxy density field (where one can form ratios of observables to separate out the 
galaxy bias). Recall that extensions to gravity act on weak lensing in multiple ways: 
the growth of the matter power spectrum alters (which can be phrased in terms of 
7), but also the light deflection law changes, involving the post-GR parameter Q 
separate from 7. This must also be included in the fit, except that many classes of 
extended gravity (such as DGP and f{R) gravity) actually have t/ = 1 on cosmo- 
logical scales. One other subtlety is that the relation of the matter power spectrum 
to the photon temperature power spectrum changes with the altered growth, mod- 
ifying the mapping between the primordial photon perturbation amplitude and 
the present mass amplitude erg- When both weak lensing and cosmic microwave 
background data are included in the constraints, the treatment of Ag and erg niust 
be made consistent. 

Figure [5] illustrates the effects of fitting the gravitational growth index 7 si- 
multaneously with determining the effective dark energy equation of state w{a) = 
Wo -f- ^0(1 — a), i.e. the expansion history. The area figure of merit in the WQ-Wa 
plane decreases by 45%, but as Hutercr & Linder 2007 pointed out the consequences 
of neglecting to fit for the growth index are worse. For an assumption mistaken by 
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Figure 2. Weak gravitational lensing probes both the growth history and expansion history 
of the universe, so failure to account for possible extended gravity effects on the growth 
overestimates the tightness of constraints on the expansion history parameters wo and Wa- 
Determination of wq and Wa simultaneously with fitting for the growth index 7 is nearly 
independent of the value of 7, here shown with typical fiducial values for GR (0.55), DGP 
gravity (0.68), and f{B) gravity (0.42). The 68% confidence level constraints on 7 are 
about 0.11 for the case shown of space-based weak gravitational lensing alone over 4000 
deg^. For the impact of other growth effects (e.g. neutrino mass, spatial curvature) on 
weak lensing see Das et al. 2011. 



A7 and weak lensing alone as a probe the derived value of Wa would be biased as 
Awa ~ 8A7. When fitting for 7 simultaneously, however, the estimation of wq and 
Wa do not depend strongly on the fiducial 7. This is a reflection of 7 being defined 
specifically to separate the expansion history influence on growth from any "beyond 
general relativity" effects (Linder 2005), making it a key method to test gravity. 
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4. Paths of Gravity 

While the model independent approach allows exploration without assuming a par- 
ticular theory, it is of interest as well to consider some specific models and their 
mapping into the post-GR parameters we have described. One can talk about three 
broad classes of extensions to gravity in terms of the physics restoring them to gen- 
eral relativity in solar system conditions, as observations require (Jain & Khoury 
2010, also see Durrer 2011, Maartens 2011, Uzan 2011 in this volume): dimensional 
reduction where below a Vainshtein radius the theory acts like GR (e.g. DGP or cas- 
cading gravity) , strong coupling where extra degrees of freedom freeze out through 
gaining a large mass in a chameleon mechanism (e.g. f{R) or scalar-tensor grav- 
ity), or screening where the extra degrees of freedom decouple and vanish through 
symmetry restoration (symmetron gravity). On cosmic scales, the first and third 
classes behave similarly, so we examine two representative cases: DGP gravity and 
f{R) gravity. 

In both of these cases, the light variable G is simply equal to unity, the GR 
value. However the growth variable V is affected. The expressions become 



Vdgp 



2 + 4r!i(a) 

3 + ^ni{a) 

3 + 3K^{k,a) ' 



(4.1) 
(4.2) 



where rim(a) — 87rG'Arpm(a)/[3i?^(a)] is the dimensionlcss matter density, H ~ d/a 
is the Hubble expansion rate, and k = k/[aM{a)] where M(a) « {2>d^ f /dR^)^^^^ 
is the effective scalar field mass. Note that DGP gravity does not have scale de- 
pendence on cosmic scales above the Vainshtein scale; gravity is scale free (e.g. the 
force is a power law with distance) on both the large scale (5-dimensional gravity) 
and small scale (4-dimensional gravity) limits and only the Vainshtein scale defined 
from the 5-d to 4-d crossover breaks this. On the other hand, f{R) gravity has both 
scale and time dependence, though tied together in a specific manner. 

Recall that the gravitational growth index 7 is related to V. For DGP gravity, 
7 = 0.68 (Lue, Scoccimarro, & Starkman 2004, Linder 2005, Linder & Cahn 2007) is 
an excellent approximation to use for calculating the matter density linear growth 
factor as a function of redshift, good to 0.2%. A mild time dependence can be 
incorporated into 7 (though this is not necessary) through Eq. (27) of Linder & 
Cahn 2007. For f{R) gravity, 7 is not generally as well approximated by a constant 
in time, and has non-negligible scale dependence at redshifts z ~ 1-3 (e.g. Tsujikawa 
et al. 2009, Motohashi, Starobinsky, & Yokoyama 2010, Appleby & Weller 2010) 
although the details depend on the specific f{R) model and parameters. 

From Eqs. (|4.ip - (j4.2p we can create phase plane diagrams of the evolution of 
the post-GR function V. This is analogous to the dark energy phase plane w-w' 
for the dark energy equation of state and its time variation, where prime denotes 
d/dlna. For dark energy, such diagrams led to clear distinction of certain physical 
classes (Caldwell & Linder 2005) as well as calibration of the compact and accurate 
parametcrisation w{a) = wq + Wa(l — a) (de Putter & Linder 2008). (Also see 
Song et al. 2010 for comparison of Q and V at a fixed time in extended gravity vs 
interacting dark energy.) 
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Figure [3] shows the results in the V-V' plane for DGP and f{R) gravity. For 
DGP gravity the equation for the phase space trajectory is 

V = -^"-(^ " (4.3) 



-i8(i-v) [y-'^ 



3V_2 

4-3V 



(4.4) 



where we used the relations 



li:„ = iw^l,n{l - ^m) ^ -3r!.„ ! ^ , (4.5) 

and r2„i is the time dependent matter density. 
For f{R) gravity, the phase space trajectory is 

V - .J^^ (4.6) 



3(1 + k2)2 

2(5^ 1)k2 

3(1 + k2)2 



-6(s-1)(V-1)(V-4/3), (4.7) 



where in the second line we parameterise M{a) — Moa"" . If we wanted to look at 
the phase space evolution with respect to inverse length k rather than time a, then 
the equation still holds, with s = 2. 

Both classes of gravity theory act as thawing cases in the nomenclature of dark 
energy phase space: they are frozen in the general relativity state (V, V) = (1,0) 
in the past, then as the Hubble parameter or Ricci scalar curvature drop from 
the cosmic expansion the theory thaws and moves away from GR. The theories 
eventually freeze to an asymptotic attractor with V" = in the future, with V = 2/3 
in the case of DGP gravity (weaker gravity) and V = 4/3 in the case of /(i?) gravity 
(stronger gravity). Note that interpreting f{R) as a scalar-tensor theory agrees with 
the expectation that gravity should strengthen, since forces carried by a scalar field 
are attractive. The clear separation of phase space territory for the two classes, 
as for thawing and freezing fields of dark energy, shows how observations could 
distinguish the nature of gravity. This also defines a science requirement that V 
should be measurable to an accuracy better than 0.1 for a 3cr distinction of gravity 
theories. 

However, there is one further point regarding f{R). While the form M{a) ^ a"^ 
is a reasonable description for the past behavior of the scalar field mass (see, e.g., 
Bertschinger & Zukin 2008, Zhao et al. 2009, Appleby & Weller 2010), in the future 
we expect M to freeze to a constant (e.g. as R itself does when the theory goes to 
the de Sitter attractor state; thanks to Stephen Appleby for pointing this out). If 
instead we parameterise the scalar field mass as M(a) = A/ia~'* + M* (accurate to 
~ 1% for at least some models), then the phase space trajectory does not asymptote 
to (4/3, 0) but rather returns to the GR limit of (1, 0) as a ^ 1 and k = k/{aM) — ^ 
in the future. 

FigureH]illustrates this behavior. As increases (dotted curve) or for wavenum- 
bers near the mass scale M (dot-dashed curve), the cosmic version of the chameleon 
mechanism begins to operate and the trajectory heads back toward GR. Also note 
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Figure 3. The phase space trajectories of DGP gravity and f{R) gravity resemble thawing 
dark energy, except here the theories move away from general relativity in the post-GR 
growth parameter plane, instead of thawing from a cosmological constant. Unlike canonical 
dark energy, the theories thaw in opposite directions: DGP moves to weaker gravity than 
GR while f{R) moves to stronger gravity. The phase space location today for each of the 
theories is shown by the squares (for Qm = 0.27). 



that because M is no longer a power law in a, different wavenumbers k do not 
simply correspond to a rescaling of a and so the single f{R) trajectory in Fig. [3] 
breaks up into varied paths for different k. The evolution along a path varies as 
well. This can be seen by the different values of V at a = 0.5 for k/AIi = 10 (blue 
square near the peak of the figure) vs k/Mi = 100 (red square near the right side). 
However, by the present (and for several e-folds to the future) this k dependence 
has vanished, with V(a — 1) shown with stars for the two cases agreeing to within 
0.5%. This holds for all k/Mi ^ 1 (note that the scalar mass Mi is likely to be of 
order the Hubble constant for observationally allowed models). Only for k/Mi « 1 
does scale dependence today enter (and the trajectory as a whole deviate). Since 
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Figure 4. In the f{R) case with a de Sitter future the scalar field mass M freezes at 
a non-zero constant, causing the theory to restore GR on cosmic scales just like with 
the chameleon mechanism on small scales. The heavy black curve recreates the M* — 0, 
vanishing mass case of Fig.[3]that freezes at V = 4/3, however the finite mass cases return 
to V = 1. While the wavenumber dependence of the post-GR parameter ai z = 1 (squares) 
is appreciable, by 2; = the fc-dependence is negligible (stars). Since the gravitational 
growth index 7 is directly related to V, this explains the scale independence of 7 at the 
present in such f{R) models. 

the gravitational growth index 7 is directly related to V, this explains the scale 
independence of 7 in such f{R) models at the present (see, e.g., Tsujikawa et al. 
2009, Appleby & Weller 2010). 



5. Conclusions 



Cosmic acceleration may be a sign that gravitation deviates from general relativity 
on large scales, pointing the way to a deeper theory of gravity and perhaps the 
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nature of spacetime and fundamental physics. Even apart from this, cosmological 
observations now have the capabihty to test gravity on scales barely probed, and 
we should certainly do so. 

In the absence of a compelling specific model, and to remain receptive to sur- 
prises in how gravity behaves, a model independent approach in terms of parame- 
terising the relations between the metric potentials and matter density and velocity 
observables has advantages. We have explored here the "2 x 2 x 2" approach pa- 
rameterising light/growth functions Q, V motivated by observable effects on photon 
and matter perturbations. These are divided into bins sensitive to both scale and 
time dependence. 

We find strong complementarity between Q and V, with each probing specific 
aspects of extensions to gravity, and both capable of being constrained by a variety 
of cosmological methods. Current surveys are making some inroads on determining 
G, but constraining V requires future large surveys such as BigBOSS, Euclid, or 
WFIRST. Crosscorrelations between different probes will be valuable as well, and 
geometric measures such as supernova and BAO distances will be essential to fit 
simultaneously the cosmic expansion history. Conversely, fitting extended gravity 
growth parameters such as 7 when using probes such as weak lensing to measure 
the equation of state is necessary to avoid bias. 

Viewing the post-GR parameters in a dynamical phase plane yields insights 
similar to its use for dark energy. Models such as DGP and f{R) gravity act as 
thawing fields, although evolving in opposite directions away from GR. Phase di- 
agrams also can illustrate the scale dependence of 7 at various epochs, and most 
importantly deliver a science requirement on distinguishing classes of gravity: fu- 
ture surveys should aim to determine V to better than 0.1. Planned next generation 
experiments such as BigBOSS, Euclid, and WFIRST can indeed potentially reach 
this level, with the caveat that expansion history should be tested as well, such as 
through supernova distances immune to gravitational modifications. 

Gravity can and should be tested on all scales, on laboratory, solar system, 
compact object, galactic, cluster, cosmological, and horizon scales. The approach 
discussed here is designed for model independence on cosmological scales, smaller 
than the horizon. On horizon scales the light/growth functions become more com- 
plicated or insufficient as other terms in the equations become important (see, e.g., 
Bcrtschinger & Zukin 2008, Hu & Sawicki 2007, Ferreira & Skordis 2010). 

Another model independent approach is to check consistency relations within 
general relativity (see, e.g., Zhang et al. 2007, Acquaviva et al. 2008, Reyes et al. 
2010). These can provide an alert to deviations, and then one must adopt specific 
models or more detailed parameterisations such as discussed here to characterise 
the physics. 
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